home *** CD-ROM | disk | FTP | other *** search
/ Info-Mac 4 / Info_Mac IV CD-ROM (Pacific HiTech Inc.)(August 1994).iso / Science / MAL ©P.f.Howden 1⁄1⁄89 / DIFFDIFFUN < prev    next >
Text File  |  1989-02-08  |  870b  |  15 lines
  1. 200 CLS:CLEAR:PRINT"DIFFDIFF FUN:":PRINT"DIFFERENTIATING A DIFFERENTIAL FUNCTION":PRINT"eg. ON LINE 230 F=X(1)^3-T^3":PRINT"FOR A 1st ORDER FUNCTION WHERE X(1)=dX/dT, X(2)=SECOND ORDER..."
  2. 201 DEFDBL B-H,O-Z
  3. 202 INPUT "MAX GIVEN FUNCTION ORDER=";N:INPUT "REQUIRED DIFFERENTIAL=";N0:N=N+N0
  4. 205 DIM X(N),Y(N):PRINT"INPUT LOWER DIFFERENTIAL VALUES:":FOR J=0 TO N:PRINT"X(";J;")=";:INPUT Y(J):NEXT J:INPUT "T=";T1:GOSUB 247
  5. 215 INPUT "INCREMENT Z (.00001?)=";Z:FOR  J=0 TO 2^N0-1:J1=J:FOR K=1 TO N0
  6. 220 IF J/2-INT(J/2)>0 THEN M=M+1:A=-1:GOSUB 250:GOTO 227
  7. 225 A=1:GOSUB 250
  8. 227 J=INT(J/2):NEXT K
  9. 230 F=X(1)^3-T^3
  10. 240 U=U+(-1)^M*F:M=0:GOSUB 247:J=J1:NEXT J
  11. 245 PRINT N0;"/DIFFERENTIAL=";U/((2*Z)^N0):GOSUB 247:U=0:PRINT"CHANGE Z?":GOTO 215
  12. 247 T=T1:FOR I=0 TO N:X(I)=Y(I):NEXT I:RETURN
  13. 250 T=T+A*Z:FOR I=0 TO N-1:X(I)=X(I)+A*Z*X(I+1):NEXT I:RETURN
  14.          MACINTOSH LISTING
  15.